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ABSTRACT 



Context. Infrared dark clouds are massive, dense clouds seen in extinction against the IR Galactic background. Many of these objects 
appear to be on the verge of star and star cluster formation. 

Aims. Our aim is to understand the physical properties of IRDCs in very early evolutionary phases. We selected the filamentary IRDC 
G351.77 - 0.51, which is remarkably IR quiet at 8yum. 

Methods. As a first step, we observed mm dust continuum emission and rotational lines of moderate and dense gas tracers to charac- 
terise different condensations along the IRDC and study the velocity field of the filament. 

Results. Our initial study confirms coherent velocity distribution along the infrared dark cloud ruling out any coincidental projection 
effects. Excellent correlation between MIR extinction, mm continuum emission and gas distribution is found. Large-scale turbulence 
and line profiles throughout the filament is indicative of a shock in this cloud. Excellent correlation between line width, and MIR 
brightness indicates turbulence driven by local star formation. 

Key words, stars; formation; stars: protostars; ISM: clouds; ISM: individual objects G35 1.77-0.51 



1. Introduction 

Although the majority of stars, especially massive stars, are ob- 
■ served to form in clusters dTesti et all 1 1 99711 1 999t iLada & LadaL 
' l2003HdeWitet al.Ll2004l) . our understanding of star formation 
• is mainly based on studies of nearby star-forming regions that 
' mostly give birth to isolated low-mass stars. On the other hand, 
little is known about the formation of clusters and of stars 
in clusters. Naturally, studies of clusters should begin by un- 
[ derstanding their youngest evolutionary phases. However, the 
. youngest (proto)clusters would still be deeply embedded within 
their natal dense clouds, thus obscured from view except at 
far IR/millimetre wavelengths. Due to the coarse angular reso- 
. lution of millimetre observations, studies of deeply embedded 
' (proto)clusters have s o far focused on the nea rest star forming 
' regions in Ophiuch tis (iMotte et al.L[T9 98: Johns tone et al.Ll2004t 



the bright Galactic background during mid-IR imaging sur- 



veys with the Infrared Space Observatory (ISO, iPerault etal 
19961) and the Mid-course Space Experiment (MSX. lEgan et al 
1998h . Their spatial coincidence with emission from molecu- 



' Andre et al., 2007"), Perseus ('johnstone et al., 2010) and Serpens 
(Testi & Sar gent. 1998; Tes ti et al, . 2000; Olmi & Testi, 2002i) . 
These regions, however, do not form massive stars and rich 
clusters and it is therefore necessary to extend these studies to 
more distant massive star-forming regions. Observational stud- 
ies of the initial conditions for the formati on of massive clusters 
have so far been v ery lim ited (e.g. Molinari et all Il998l I2002L 
iBeuther & Sc"hiika. I200I but will soon find a new renaissance 
with the high angular resolution and sensitivity offered by the 
Atacama Large Millimeter/sub millimeter Array (ALMA). 

Infrared dark clouds (IRDCs) are ideal candidates for study- 
ing the earliest phases of the formation of clusters. They were 
discovered as dark patches in the sky at IR wavelengths against 
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lar lines and dust indicates that they consist of dense molecu- 
lar gas and dust of a few 10^ cm""' density (Carey et al., 199^ 
200 0D. Ternperatur es in these clouds are low enough (10-20 K, 
Pillai et al.L l2006al) that they do not radiate significantly even at 
mi d-IR wavelength s : of th e 190 dust cores in IRDCs surveyed 
by 'Rath borne et al.l (1201 Ol) only 93 have detectable emission at 
24jum. 

The forma tion mechanism of IR DC is still a matter of debate 
(lMyers','2009^). Re cent r esults (e.g., i Jimenez-Serra et al.Ll2010l; 
Hernandez & TanLl201 iT) may support models of filament forma- 
tion from converging flows. However, there is still no solid con- 
clusion about their origin. Tan (2007) suggests that IRDC may 
contain the sites of future massive star and star cluster formation, 
since their densities and mass surface densities are similar to re- 
gions known to be undergoing such formation activity. Indeed, 
emerging evidence suggests that the earliest phases of the for- 
mation of hi gh-mass star and cluster are occurring wit hin these 
clouds (e.g.. iPiUai et"an.l2006btlRathborne et alll2006l) . Studies 
of IRDCs have often focused on sources that already have in- 
tense star formation activit y and are therefore in late protostellar 
evolutionary phases (e.g., IPillai et al.L l200 6b; Rathborn e et all 
2OO81) . However, to address questions related to the formation of 
clusters, and of the single stars in them, studies of earlier evolu- 
tionary phases with no or modest star formation are needed. 

The IRDC under investigation in this paper is at Galactic 
coordinates (/>]=(35r.77,-00°.51), and was identified by 
ISimon et all (l2006l) in the 8jum data from the MSX satelUte as 
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G351.77 - 0.51 (hereafter G351.77). G351.77 appears as a fil- 
ament of dust, which is seen in absorption against the mid-IR 
backg round of our Galaxy (see Fig. [TJ up to lAjim (Carey et al., 
2009 ) and in emi ssion in the millimetre range (iFaundez et all 
20041) . Following iMotte et al.1 (HoOTj), we used the 8/zm MSX 
image as a starting point to infer the stellar activity of the fila- 
ment. G35 1 .77 is remarkably quiet at this wavelength: emission 
is detected towards five positions in the filament, three of them 
being close to the bright IR source IRAS 17233-3606. Compact 
emission is detected in more sensitive Spitzer images towards 
other positions in the filament (see Sect. |3]l. IRAS 17233-3606 
lies near one end of the filament. This source shows sev- 
eral signs of active massive s t ar formation: differen t authors 
(ICaswell et all [T980l iFixetaLl ITosl iMeiitenl [1991 detected 
very intense H2O, OH, and CH3OH masers, while i Leurini et all 
(12008. 2009) report the discovery of bipolar outflows originat- 
ing in the vicinity of IRAS 17233-3606 and a molecular spec- 
trum typical of hot molecular cores near massive young stellar 
objects. If IRAS 17233-3606 is part of the molecular environ- 
ment of G351.77, the distance of G351.77 would be relatively 
clos e, since studies of I RAS 17233-3606 locate it at D < 1 kpc 
(see iLeurini et all |201 1 . and reference therein). The near dis- 
tance would imply that high linear resolution observations could 
be performed even with single-dish radio telescopes whose typi- 
cal angular resolutions are not better than 10", coiTesponding to 
0.05 pc at 1 kpc. In addition, given its low declination, G351.77 
would be an ideal candidate for high angular resolution obser- 
vations with ALMA to study the earliest stages of formation of 
massive stars and stellar cluster, and to investigate the content of 
low- and intermediate-mass star formation in IRDCs. 



Table 1. Summary of the parameters of the APEX and MOPRA 
line observations. 



Setup 


Tuning Frequency 


Telescope 


Beam 


Av 




(MHz) 




(") 


(km s-') 


N2H+ (1-0) 


93173.5 


MOPRA 


35'.'0 


0.1 


'^CO-C'^O (2-1) 


219980.0 


APEX 


28'.'4 


0.17 


C"0 (2-1) 


224700.0 


APEX 


27'.' 8 


0.16 



Table 2. Positions observed with APEX 



Position 


R.A. [J2000] 


Dec. [J2000] 


Setup 






Clump- P' 


17:26:42.30 


-36:09:18.23 


"CO-C'O, 


c" 





Clump-2 


17:26:38.79 


-36:08:05.53 


"co-c'*'o, 


c" 





Clump-3 


17:26:47.33 


-36:12:14.37 


i3C0-C'*'0, 


c" 





Clump-4 


17:26:40.30 


-36:08:29.63 


'^CO-C' 


^0 




Clump-5 


17:26:25.25 


-36:04:57.03 


'3CO-C'**0, 


c" 





Clump-6 


17:26:23.25 


-36:04:32.74 


13CO-C' 


^0 




Clump-7 


17:26:34.78 


-36:06:34.22 


'3C0-C'**0, 


c" 





Clump-8 


17:26:28.26 


-36:05:21.34 


'^CO-C 






Clump-9'' 


17:25:31.96 


-36:22:15.70 


'3C0-C'*'0, 


c" 





Clump- 10^ 


17:26:06.66 


-36:22:09.63 


'3C0-C'**0, 


c" 





Clump- 12'' 


17:25:31.97 


-36:21:21.05 


'^CO-C 


^0 




G351.62 


17:26:34.80 


-36:17:41.00 


'^CO-C' 


^0 




" IRAS 17233- 


-3606 











Clump-9 and 12 are offset from IRAS 17221-3619 
IRAS 17227-3619 



2. Observations 

2.1. APE)^ observations 

The Large Apex BOlomoter CAmera (LABOCA.'S iringo et al.L 
[2009) was used to observe the continuum emission at 870//m 
towards the dust filament G351.77. The FWHM beam size of 
APEX at 870//mis ~ 19'.'2. The data are part of th e ATLASGAL 
survey of the Galactic plane (ISchuUer et al.L[2009h . Observations 
were done with fast-scanning, linear, on-the-fly maps in total- 
power mode (i.e. without chopping secondary). Details on the 
data reduction are given in iSchulleret al.l (l2009h . In particular, 
the subtraction of correlated skynoise results in a filtering of 
uniform extended emission on scales larger than ~2.5'; compact 
or filamentary sources should not be strongly affected, though. 
For a d etailed discussion of fi ltering of extended structures, we 
refer to iBelloche et all (l201lh . The calibration of the LABOCA 
data relies on opacity correction based on skydips and on ob- 
servations of planets as primary calibrators and of well-known, 
bright, compact sources as secondary calibrators. As a result, the 
photometric uncertainty is of the order of 10%. 

Observations in ^^CO, C"*0 and C'^O (2 - 1) towards se- 
lected positions in G351.77 were performed during the science 
verification programme (project number 081.F-9805(A)) of the 
APEX-1 receiver on APEX in 2008 June (19 and 25) and August 
(21, 30 and 31). The observations in June were performed under 
mediocre weather conditions (with 2-3 mm precipitable water 



'Based on observations collected at the European Organisation 
for Astronomical Research in the Southern Hemisphere, Chile, under 
programme ID 081.F-9805(A). APEX is a collaboration between 
the Max-Planck-Institut fiir Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 



vapour), yielding system temperatures of 240-400 K. The obser- 
vations in August were, on the other hand, performed under ex- 
cellent atmospheric transmissions (down to 0.4 mm precipitable 
water vapour), yielding system temperatures of 160-210 K. The 
observed transitions and basic observational parameters are sum- 
marised in Table [Tj 

For t he observations, we used the APEX facility EFT spec- 
trometer dKIein et all [20061) . The '^CO and C"*0 (2 - 1) tran- 
sitions were covered in one frequency setup, with the APEX-1 
receiver tuned to 219.98 GHz in the lower side band. A sec- 
ond frequency setup was used to cover the C'^0 (2 - 1) line, 
with the APEX-1 receiver tuned to 224.7 GHz in the lower 
side band. For both setups, we performed long integration ob- 
servations on selected positions (Table |2]i. For the '^CO-C'^O 
setup, we observed an additional strip of five points ((20", -20"), 
(10",-10"), (0",0"), (-10",10"), and (-20",20")), centred on 
the position G351.62 of Table|2] For both frequency setups, data 
were taken in the position-switching mode, with a reference po- 
sition (600", 0") in horizontal coordinates from the targets. 

We used a main-beam efficiency of 0.7^ to convert antenna 
temperatures into main-beam temperatures. A detailed descrig- 
tionof APEX and of its performance is given by iGiisten et al] 
([2OO6I) . 



^http://www.apex-telescope.org/telescope/efficiency/index.php 
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Fig. 1. Image of the S/im (left) and 24/im (right) emission from Spitzer of the molecular cloud G35 1 .77 - 0.5 1 . The numbers label 
the 870/jm dust continuum clumps found by the CLUMPFIND procedure and discussed in the text (see TableO. 



2.2. /WOPR/Q observations 

The observations were carried out with the ATNF Mopra tele- 
scope on 2007 June 30 under clear sky conditions. Two adjacent 
fields were observed in 220 x 220 arcsec^ on-the-fly maps for 30 
minutes each. The HEMT receiver was tuned to 92.2 GHz and 
the MOPS spectrometer was used in narrow band mode to cover 
the 8 GHz bandwidth of the receiver with about 0. 1 km/s velocity 
resolution. The beam size at these frequencies is about 35 arc- 
sec. The system temperature was 175 K. Pointing was checked 
and corrected with line pointings on the SiO maser source AH 
Sco. 

Initial processing of the data was done using the ATNF 
tools livedata and gridzilla and included the subtraction of the 
off signal, the calibration to T* units, the baseline subtraction 
and the gridding. To convert antenna to main beam brightness 
teruperature, we corrected the data for an antenna efficiency of 
0.5d Detai ls on antenna efficiency in the 3mm band is given by 
iLadd et all (|2005). The data was then exported to MIRIAD and 
GILDAS for further processing. 



3. IR emission 



Figures [T(a)| and [T(b)| show G351.77 at 8 and 24yum, respectively, 
from the GLIMSPE and MIPSGAL sur veys o f the Galactic 
plane (iBeniamin et al.L 120031: ICarev et al.L l2009l) . The filament 



'The MOPRA telescope is part of the Australia Telescope which 
is funded by the Commonwealth of Australia for operation as a 
National Facility managed by CSIRO. The University of New South 
Wales Digital Filter Bank used for the observations with the MOPRA 
Telescope was provided with support from the Australian Research 
Council 

''http://www.narrabri.atnf.csiro.au/mopra/obsinfo.html 



is dark at both wavelengths except for bright emission asso- 
ciated with IRAS 17233-3606 (dust condensation 1, hereafter 
clump 1) and with the dust condensation 3 (clump 3) to the 
south of IRAS 17233-3606 (see Table O. Additionally, weak 
emission, which is still compact in most of the cases, is de- 
tected at several positions along the filament at both wave- 
lengths. Compact emission at 24jum is an indicator of active star 
formation, since this emission traces warm dust heated as ma- 
terial accretes from a core onto a central protostar To search 
for evidence of outflow activity from young stellar objects, 
we also checked the 4.5//m map of the region made with the 
Spitzer InfraRed Array Camera (IRAC). Extended green emis- 
sion (4.5yum) is clearly associated only with the outflows orig- 
inating in IRAS 17233-3606 (Leurini et al., 2009). A jet-like 
feature extends from IRAS 17233-3606 into the dark filament. 
More compact emission at 4.5/im is also detected towards three 
additional positions in G351.77 (clumps 2, 3, and 5). 

A larger map of the region in the IRAC 8;um image reveals a 
complex picture: extended emission along a semicircular struc- 
ture is detected opposite to G351.77. Moreover, the map is dom- 
inated by the infrared source IRAS 17221-36 19, which is asso- 
ciated with bubble CS84 of the catalogue of IChurchwell et aP 
(1200 7). and belongs to a second circular structure seen in emis- 
sion in the mid-IR (Fig. [3]). Several dark patches are seen at 8yum 
and 24fj.m, two of them being the IRDCs G351.62 - 0.58 and 
G351.7 7 - 0.60 (see Fig. [3] and Table |2li identified in the cata- 
logue of lSimonetaP ([2006). 

4. Continuum emission at 870jum 

The continuum emission at 870/^m probes cold dust. Therefore, 
it can be used to locate dense compact clumps, potential sites 
of star and cluster formation. Given the linear resolution of our 
LABOCA map (19'.'2 ~ 0.1 pc at 1 kpc), our data are sensi- 
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Fig. 2. LABOCA map of the SlOfim dust continuum emission 
from the IRDC G35 1.77-0.51. The white contours are from 5cr 
(0.425 Jy/beam) to 2 Jy/beam in steps of 5cr; the black contours 
are from 2 Jy/beam to 5 Jy/beam in steps of lOcr. The num- 
bers label the clumps found by the CLUMPFIND procedure and 
discussed in the text (see TableO. 



five to dust clumps possibly associated with cluster format ion 
(IWimams et al.L 1200(1 iTbsti et al.Ll2000tlblmi & TestiLlIOOl . 

Our LABOCA map of G35 1 .77 reveals extended continuum 
dust emission that matches the morphology of the mid-IR ex- 
tinction extremely well (Figs. |2] and [S]). The IRDC G351.77 - 
0.60 coincides with the 6cr LABOCA emission on the east 
of IRAS 17233-3606 (Fig. |3]l. On larger scales, the emission 
extends over a curved structure (Fig. O, which includes the 
IRDC G351.62 - 0.58 (see previous section) associated with 
our position G351.62, and the sources IRAS 17227-3619 and 
IRAS 17221-3619 (see Table |2]i. A second, intense peak of 
continuum emission is found opposite to IRAS 17233-3606, 
on its west. This emission is, however, associated with the 
IR source IRAS 17220-36 09, whose molecular lines have a 
velocity of -98.5 km s ' dBronfman et a l.. 1996), which im- 
plies a kinematic distance of 10 kpc dFaundez et al., 200j). 
Since IRAS 17233-3606 has a velocity of -3.4 km s^', the 
two sources are unlikely to be physically related. We conclude 
that IRAS 17220-3609 is projected close to the Hne of sight of 
IRAS 17233-3606, but it is much further away from the Sun, so 
we did not include IRAS 17220-3609 in our following analysis. 

No continuum emission at 870;um is associated with the 
semicircular structure opposite to G351.77 seen at 8//m (Fig. [3]), 
hence probably due to low density gas not sampled by our ob- 
servations. Since our molecular line follow-up observations were 
only performed toward peaks of emission at 870;um, we do not 
have any information on the velocity field of this structure and 
we cannot determine whether it is physically associated with 
G351.77ornot. 



Table 3. Results of the CLUMPFIND calculations 



Clump 


Sy 


Fy 






/lb 

9 






[Jy] 


[Jy/beam] 


["] 


["] 


["] 








G351.77 -0.51 






1 


159.4 


54.2 


43.6 


38.5 


41.1 


8, 24yum 


2 


29.7 


6.2 


31.6 


46.7 


39.3 


8, 24yum 


3 


24.9 


5.2 


31.5 


49.6 


40.7 


8, 24yum'' 


4 


16.1 


4.0 


46.9 


25.2 


36.5 




5 


9.3 


3.1 


18.8 


20.9 


19.8 


8, 24yum 


*- 




6.1 


2.3 


21.5 


19.3 


20.4 


24;um 


7 


8.7 


2.1 


33.0 


41.5 


37.3 




8 


4.6 


1.5 


17.5 


20.5 


19.1 


24^m 


11 


3.2 


1.1 


20.6 


22.7 


21.7 


8, 24/zm 






other positions 








9 


4.8 


1.5 


28.9 


21.1 


25.1 


8, 24;um 


10 


2.7 


1.2 


23.1 


12.4 


18.2 


8, 24;um 


12 


5.0 


1.0 


46.9 


21.5 


34.8 


8, 24;um 



" deconvolved FWHM, and FWHM,. 

^ deconvolved geometric mean of FWHM^ and FWHMy 

^ association with IR emission 

the IR source is offset from the mm source 

position contaminated by the emission from IRAS 17233-3606 



4.1. Identification of tlie clumps 

To identify clumps in the millimetre continuum emission and 
define their properties, we used a two-dimensional variation 
of the clump- finding algorithm CLUMPFIND developed by 
IWilliams et al.l (1199 4). The algorithm works by eff'ectively con- 
touring the data at a multiple of the rms noise of the map, 
then searching for peaks of emission to locate the clumps, and 
finally following the clump profile down to lower intensities. 
The advantage of this method is that it does not a priori as- 
sume any shape for the clumps. However, since the contour- 
ing levels are chosen by hand, the algorithm does not take 
any emis sion into account b elow the user-selected threshold. 
Recently, iPineda et al.l (12009 ') have outlined some weaknesses 
of the method, which may lead to erroneous mass functions 
for crowded regions. Nevertheless, the CLUMPFIND algorithm 
may still be useful for studying the structure of one cloud as in 
our case. We set the threshold level to 5cr. 

The procedure calculates the peak position, the full width 
at half maximum not corrected for beam size for the x-axis, 
FWHM,, and for the y-axis, FWHMy, and the total flux density 
integrated within the clump boundary within the threshold level. 

We identified twelve clumps in the 870yL(m continuum emis- 
sion over the region shown in Fig. |3] The centroid coordinates 
dehvered by CLUMPFIND are given in Table |2] Other results 
of the CLUMPFIND analysis are given in Table [3] where we 
also indicate whether clumps are associated with continuum 
emission at 8 and 24yum, either point-like or diffuse. We inter- 
pret the detection of a compact 24yum source in all clumps of 
G35 1 .77 except clumps 4 and 7 as a confirmation of the validity 
of the method used for identifying clumps in the region. While 
clump 7 is clearly detected as a secondary peak of emission at 
870//m, clump 4 does not look like a clumpy structure at vi- 
sual inspection, and could be an example of 'pathological' fea- 
tures created by CLUMPFIND in case of significant emissio n 
found between prominent clumps (see iGoodman et al.L l2009l) . 
Additionally, we note that clump 1 is associated with the ac- 
tive massive star-forming region IRAS 17233-3606 discussed 
above, and that clumps 9 and 12 are associated with the in- 
frared source IRAS 17221-3619, which hosts an Hii region (e.g.. 
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Fig. 3. Left: Continuum emission at 870/vm of the environment around G351.77 - 0.51. The numbers label the clumps found by 
the CLUMPFIND algorithm outside the region plotted in Fig.[T] The white cross marks the position G351.62, the stars the IRAS 
sources found in the region. Right: Image of the emission at 8jum of the molecular environment around G35 1 .77 - 0.5 1 . The white 
(and black for IRAS 17221-3619) contours show the emission at 870/im from 0.255 (3 cr) to 4 Jy/beam in step of 10 cr. In both 
panels, the white box outlines the region plotted in Fig.[T] 



Table 4. Physical parameters of the clumps derived from the LABOCA emission. 



Clump 


MlOK 


M25K 


M35K 




-^H2.25K 


^H2,35K 


"Ht.IOK 


"H2,25K 


"H2,35K 




[Mo] 


[Mo] 


[Mo] 


[10^2 cm-2] 


[10^2 cm--] [10-2 cm-2] 


[10^ cm-3] 


[10^ cm-3] 


[10^ cm-3] 












G351.77 -0.51 










1 


2964 


664 


428 


460 


103 


66 


108 


24 


16 


2 


552 


124 


80 


52 


12 


8 


23 


5 


3 


3 


463 


104 


67 


44 


10 


6 


17 


4 


3 


4 


300 


67 


43 


34 


8 


5 


16 


3 


2 


5 


173 


39 


25 


27 


6 


4 


56 


13 


8 


6 


113 


25 


16 


19 


4 


3 


34 


7 


5 


7 


161 


36 


23 


18 


4 


3 


8 


2 


1 


8 


85 


19 


12 


13 


3 


2 


31 


7 


4 


11 


60 


13 


9 


10 


2 

other positions 


1 


15 


3 


2 


9 


89 


20 


13 


13 


3 


2 


14 


3 


2 


10 


51 


11 


7 


10 


2 


1 


21 


5 


3 


12 


94 


21 


14 


9 


2 


1 


6 


1 


1 



iMartm-H ernandez et al.', '2003*) surrounding an IR star cluster 
candidate (Borissova et al.. .2006.) associated with bubble CS84. 

4.2. Mass estimate 

The masses of the clumps given in Table|4]are derived by assum- 
ing that the continuum emission at 870yL(m is optically thin: 



KyBy(Td) 



(1) 



where Fy is the integrated flux, d the distance of the source, k the 
dust opacity, and BiTj) is the Planck function for a black body 
of dust temperature T^, all for a frequency of ~ 340 GHz. For 
the dust opacity, we adopted a value of 0.0182 cm^ g"' (from 
Table 1, Col. 5 of Ossenkopf & Henning, 1994). This agrees 
with recent measuremen ts of dust opacities at 850 and 450yum 
from lShirlev et al.l (1201 ih in the low-mass Class core B335. 



From the line observations performed towards the dust fil- 
ament of G351.77, we confirm that IRAS 17233-3606 belongs 
to the same molecular environment of the infrared dark cloud 
G351.77 (see Sect.|5]i. We assume a distance of 1 kpc in the fol- 
lowing calculations. Our mass estimates would be a factor of 2 
smaller for a distance of 700 pc (see Sect. [1]), corresponding to 
the s mallest distance o f IRAS 17233-3606 reported in the liter- 
ature (iJVIiettinen et al.L 12006). We also computed beam-averaged 
H2 volume densities assuming that the clumps have spherical 
symmetry and a mean molec ular weight per hy drogen molecule, 
jUHt, of 2.8 (see discussion in lKauffmann et al] l2008). From pre- 
vious line observations, w e know that c l ump 1 (IRAS 17233- 
3606) harbours a hot core dLeurini et al.L l2008l) . while our cur- 
rent data towards the other clumps do not reveal any rich molec- 
ular spectrum indicative of hot core activity, i.e., an embedded 
central heating source. Ammonia (1,1) and (2,2) observations 
with the Parkes telescope (Wienen et al. in prep.) of three posi- 
tions along the filament (clump 1, 2, and 3) infer temperatures 
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of <20 K on an angular size of ~ 1'. Except for clump 1, which 
hosts the hot core associated with IRAS 17233-3606, we believe 
that temperatures of 10-25 K are appropriate for the dust in the 
clumps in G351.77 where only point-like emission is detected 
at 24jum. We computed masses, column densities, and densities 
of all clumps also for Td = 35 K for consistency with Sect. 15.21 
Results can be easily scaled to higher temperatures: compared 
to the values given for Td - 35 K, estimates are a factor ~ 1 .5 
smaller for Td - 50 K, likely to be an upper limit to the temper- 
ature of the dust for IRAS 17233-3606 at the resolution of the 
data. For T = 25 K, the lowest mass derived from our analysis 
for the clumps in G3512.77 is 13 Mq (clump 1 1), while all other 
clumps have masses of several tens of solar masses. As a result, 
G351.77 has the potential of forming intermediate to high-mass 
stars or clusters of low-mass stars in several positions along the 
filament. 



5. Emission from CO isotopologues 

5.1. Velocity Held 

To derive the velocity of each clump, we used the optically thin 
C"0 (2 - 1) line (see Sect. EIll. The C'**0 (2 - 1) transition 
was used for the positions that are not observed in C'^O (see 
Table|2]). Both transitions were detected towards all positions to- 
wards which they were observed. To derive the parameters of 
the C"0 line, we used the HFS method in CLASS. A Gaussian 
fit was applied to the C'**0 and '^CO (2 - 1) Unes, although 
their spectra show more complex velocity structures. The '^CO 
line often shows evidence of non-Gaussian wings and multiple 
peaks. These are probably due to self absorption, as the optical 
thinner C'^O and C'^O transitions peak at the velocity of the dip 
in the '■'CO spectrum. Negative features at other velocities are 
most likely due to emission in the reference positions used dur- 
ing the observations. However, in the clumps that were observed 
in all three isotopologues, no absorption feature is detected in the 
C'^O spectrum. Therefore we believe that our results (based on 
the analysis of the C'^O spectra, see Sect. 16.1b are not affected 
by this problem. The spectra of all clumps are shown in Figs.|4] 
and|5] while the results of the fit are listed in Table|5] 

From the analysis of the velocity field, we found two dif- 
ferent molecular components along the ring of dust detected at 
870 yum, and all clumps along the infrared dark cloud G35 1 .77 
are at velocities -3 km s"', including clump 3, immediately 
south of IRAS 17233-3606 (clump 1). We confirm the 5cr 
dust continuum detection at the position G351.62 (q'2oo() - 
17''26"'34'.80, 52000 = -36°17'41'.'00) to be due to dense mate- 
rial, since both '^CO (2 - 1) and C'**0 (2 - 1) are detected at that 
position, with a velocity similar to that of G351.77. Clumps 9 
and 12 (see Fig.[3]l have velocities around -22 km s"' . Moreover, 
at three positions (clump 3, clump 10, and G351.62) two ve- 
locity components are detected in emission, in '^CO(2 - 1) at 
clumps 3 and 10, and in CO (2 - 1) and C'**0 (2 - 1) at all 
five points across the position G35 1 .62. The two velocity com- 
ponents (clump 5 -3.2 and -24.5 km s"'; clump 10: ~ -4.0 
and -21.5 km s"'; G351.62: -2.6 and -21.9 km s"') belong 
to the same molecular environment of the infrared dark cloud 
G351.77 and to that of clumps 9 and 12. On clump 4, the feature 
at Vlsr - 3.2 km s"' shows non-Gaussian blue-shifted wings; 
on clump 10, the feature at -21.5 km s ' also has non Gaussian 
red-shifted wings. 



^velocities of clumps 3 and 10 are here derived from the '^CO (2-1) 
line since the second component is not detected in 



5.2. Column density estimate 

From the integrated intensity of the C'^O (2 - 1) line, or al- 
ternatively from the C'^O (2 - 1) transition for those sources 
not observed in C'^O, we can derive the beam-averaged col- 
umn density of CO, assuming a given abundance of CO rel- 
ative to its rarer isotopolog ues (Xieo/isQ ~ 560, XuQ/nQ ~ 
1800. IWilson & RoodL Il994l) . and a given excitation tempera- 
ture. Although the HFS method in CLASS infers the opacities 
for each components of the hyperfine structure of a transition, 
the results are very uncertain for overlapping components. This 
is the case for the C'^O (2-1) line, where the largest separation 
in velocity between the diff'erent components is ~ 2.4 km s"' 
and lines have typical widths of ~ 2 km s '. 

We therefore computed the C'^O column density under the 
assumption tha t the (2-1) line i s optically thin and that the 
gas is in LTE dCaselU et al.L |2002|) . iHofner et al.l (|2000|) found 
rotational temperatures for C'^O between 16 and 41 K in a sam- 
ple of ultra-compact HII regions. Therefore, we computed col- 
umn densities for Tex = 10, 25 and 35 K. The partition function 
QiTex) was estimated as 2(^ex) - aT^, where a and /3 are the 
best-fit parameters from a fit to the p artition function obtained 
from CDMS catalogue jMiilleret aU i2001) at diff'erent excita- 
tion temperatures in the range 10-500 K. The values of Q{Tg^) 
at 10, 25, and 35 K are 3.8, 9.4, and 13.2, respectively for C'^O, 
3.9, 9.6 and 13.5 for C'^O. The results are given in Tables |6] 
and|2l The uncertainties are statistical errors based on the fit (see 
Table |5]i and do not take the uncertainties on the excitation tem- 
perature into account. We refer the reader to Sect. 16. II where the 
effect of an uncertain excitation temperature on column density 
is considered when estimating molecular abundances. Moreover, 
the column densities derived for the two isotopologues of CO 
can be lower limits to the real values in case the lines are opti- 
cally thick. 

For the C"^0(2 - 1) line we can verify if the assump- 
tion of optically thin emission is correct for those clumps that 
have been observed in both isotopologues. The expected ratio 
^18,17 ijetYyeen the integrated intensities of the two transitions 
should be similar to th e canonical value of 3.65 (Penzias, 198T|; 
IWilson & Rood 1 19941) for the abundance of C'^O relative to 
C'^O. We found ^'^"'^ > 3 for clumps 7 and 10, while clumps 
1, 2, 3, and 9 have values of ~ 2.5, which implies that the C'^O 
(2-1) transition is partially optically thick. In addition, clump 5 
has ~ 5.3, which could imply that the two isotopologues 

may have different excitations or extent of emission. We there- 
fore conclude that the column densities derived from the C'^O 
transition are likely lower limits to the real values. 



6. Emission from N2H+ 

The lower excitation lines of NaH^ have been shown to be 
excellent tracers of cold and dense gas (iTafalla et alJ l2004 
iBergin et alj|2002l) . With the goal of deriving the properties of 
the dense gas, we mapped the infrared dark patch of G35 1 .77 in 
N2H^ (1-0). As expected, dust and N2H^ integrated emission 
correlate very well (see Fig. |6|l. Though the resolution of the 
MOPRA data is approximately two times lower than that of our 
dust continuum data, at least three bright cores can be identified 
along the filament. These coincide with the dust clumps labelled 
1, 2, and 5. Clumps 4, 6, 7, and 8 are also associated with promi- 
nent N2H^ emission, although not clearly resolved into cores. 
Clump 1 1 is associated with a ridge of low-level N2H^ emis- 
sion, consistent with the low dust mass of this clump. 
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Fig. 4. In grey, spectra of '^CO (2-1) of clumps 1, 2, 3, 4, 5 and 6. The solid and dashed spectra are C"^0 (2-1), and C'^'O (2-1) 
respectively. 
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Fig. 5. In grey, spectra of '"'CO (2-1) of clumps 7, 8, 9, 10, 12 and position G351.62. The solid and dashed spectra are C'**© (2-1), 
and C' ''O (2-1) respectively. 



The interaction of the molecular electric field gradient 
and the electric quadrupole moments of the nitrogen nuclei 
(1=1) leads to hyperfine structure and a splitting of the N2H^ 
(1-0) transition into seven components. The large line widths 
often attributed to the high level of turbulence in high-mass 
star-forming regions, lead to blending and result in three 
blended groups (one main and two satellites). As in the case 
of C'^O, the CLASS HFS method has been used to fit the 
line profile. Thus, we can accurately determine the line optical 
depth, excitation temperature, line width, and LSR velocity. 
These parameters were determined for the four clumps averaged 
over the MOPRA beam of 35". Subsequently, we determined 



the line column density following ' Caselli et al.l (l2002l) . The 
spectra of the three clumps are shown in Fig. |7] while the 
line parameters are tabulated in Table [8] The heavy blending 
of hyperfine components means less stringent constraints on 
the optical depth and excitation temperature. In particular, the 
optical depth for clump 1 derived from the HFS method in 
CLASS is unexpectedly low (0.2) with > 50% eiTor. The fit to 
the hyperfine structure is as good for a range in optical depth, 
0.2 < Ttot < L Hence, in Table |8] we give the range in column 
densities, abundance and abundance ratio bracketing Ttot - 0.2 
and Ttot = 1 ■ Fits for Ttot > 1 were poorer and discarded. 
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Table 5. Line parameters of the C'^O and C'**0(2 - 1) transitions. The symbol - means that the line was not observed. 



Clump 


Vlsr 


Av 






Vlsr 


Av 


7'~^ 7- 

J T^bSy 




[km s-'] 


[km s"' 


] 


[Kkms-'] 


[km s-'] 


[km s-'] 


[Kkms-i] 






C"0 








C'«0 




1" 


-2.51 ±0.04 


5.0 ±0, 


1 


32.1 ±0.5 


-2.65 ± 0.02 


5.68 ± 0.04 


84.0 ± 0.5 


2 


-2.21 ±0.02 


2.4 ± 0, 


1 


8.4 ±0.1 


-2.28 ± 0.02 


3.33 ±0.06 


21.3 ±0.3 


3 


-2.95 ± 0.06 


1.4 ±0, 


,2 


4.2 ±0.2 


-3.16 ±0.02 


2.17 ±0.05 


10.1 ±0.2 


4 


- 


- 




- 


-3.16 ±0.02 


3.66 ± 0.06 


20.8 ±0.3 


5 


-2.6 ± 0.1 


2.2 ± 0, 


.3 


2.3 ± 0.2 


-2.66 ± 0.03 


2.82 ± 0.07 


11.7 ± 0.3 


6= 










-2.5 ±0.1 


4.0 ±0.3 


6.9 ±0.4 


7 


-3.35 ± 0.03 


1.9 ±0, 


,2 


3.52 ± 0.07 


-3.19 ±0.05 


2.6 ±0.1 


12.1 ±0.4 


8" 










-3.39 ±0.05 


2.5 ±0.1 


6.8 ±0.3 


9 


-21.80 ±0.02 


1.94 ±0, 


.09 


6.33 ± 0.08 


-21.70 ±0.02 


2.43 ± 0.05 


14.9 ±0.3 


10 


-22.80 ± 0.07 


2.1 ±0, 


,2 


1.39 ±0.08 


-22.67 ± 0.08 


1.9 ±0.2 


4.5 ±0.3 


12^ 










-22.27 ± 0.04 


3.1 ±0.1 


11.1 ±0.3 


G351.62f 










-2.58 ±0.03 


1.66 ±0.07 


4.2 ±0.2 












-22.0 ±0.1 


1.3 ±0.2 


0.8 ±0.1 



Uncertainties are the errors in the Gaussian and/or HFS fit, and do not include calibration uncertainties. 

" the C"0(2 - 1) integrated intensity has been calculated from a Gaussian fit. 

the C"*0(2 - 1) line has a non Gaussian red- and blue-shifted emission. 

the C'**0(2 - 1) line has a double peak profile. The fit was performed with one component 

the C"*0(2 - 1) line has a non Gaussian red-shifted emission. 
^ the C"*0(2 - 1) line has a non Gaussian blue-shifted emission. 

' the C'*0(2 - 1) line has two velocity components (as well as the ''CO (2 - 1) transition); the one at Vlsr = -2.6 has non Gaussian blue-shifted 
emission. 



Table 6. Column densities of C'^O. 



6.1. CO-N9H+ enhancement 



Clump 




(10'"* cm-2) 








Tex = 10 K 


Tex = 25 K 


ex 


35 K 


1 


38.8 ±0.6 


35.6 ±0.6 


41.3 = 


tO.6 


2 


10.2 ±0.1 


9.3 ±0.1 


10.8 = 


tO.l 


3 


5.1 ±0.2 


4.7 ±0.2 


5.4 = 


tO.3 


5 


2.8 ±0.2 


2.6 ±0.2 


3.0 = 


tO.3 


7 


4.25 ± 0.08 


3.91 ± 0.08 


4.5 = 


tO.l 


9 


7.7 ±0.1 


7.03 ± 0.09 


8.1 = 


tO.l 


10 


1.7 ±0.1 


1.54 ±0.09 


1.8 = 


tO.l 


7. Column densities of C 








Clump 




A'ciso 
(10'5 cm-2) 








T,, = 10 K 


Tex = 25 K 


Tex = 


35 K 


1 


104.0 ± 0.6 


97.5 ± 0.6 


113.4 


±0.7 


2 


26.4 ± 0.4 


24.7 ± 0.3 


29.6 


±0.4 


3 


12.5 ±0.2 


11.7 ±0.2 


13.6 


±0.3 


4 


25.6 ±0.3 


24.1 ±0.3 


28.1 


±0.4 


5 


14.5 ±0.5 


13.6 ±0.5 


15.8 


±0.5 


6 


8.5 ±0.5 


8.0 ±0.5 


9.3 


±0.5 


7 


15.0 ±0.5 


14.0 ±0.5 


16.3 


±0.5 


8 


8.4 ±0.4 


7.9 ±0.3 


9.2 


±0.4 


9 


18.4 ±0.4 


17.3 ±0.3 


20.1 


±0.5 


10 


5.6 ±0.4 


5.2 ± 0.3 


6.1 


±0.5 


12 


13.7 ±0.4 


12.9 ± 0.3 


15.0 


±0.5 


G351.62 


5.2 ±0.2 


4.9 ± 0.2 


5.7 


±0.3 




1.0 ±0.1 


0.9 ±0.1 


1.1 


±0.1 



For reasons not yet fully understood, N2H^ (along with am- 
monia), unl ike CO, is resistant to depletion in very dense and 
cold cores dOberg et al.L l2005h . However, internal heating re- 
introduces CO to the gas phase. Since CO is one of the main 
destroyers of N2H^ in the gas phase (e.g., Aikawa et al., 200 1|, 
the abundance of N2H^ decreases when protostellar/cluster for- 
mation heats up the envelope to sufficiently high temperatures, 
releasing the frozen out CO. The ratio of the observed CO- 
to-N2H^ abundance should thus reflect this chemistry. Hence, 
this differential depletion of the two species might therefore be 
used as a chemical clock to assess the evolutionary status of the 
clumps along the filament. 

After smoothing the LABOCA data to the same resolution 
as the N2H^ and C'^O data, we estimated the molecular abun- 
dances of the two species and the relative abundance of C'^O 
to N2H^. Results are tabulated in Table [8] For this, we used 
the C'^O column density assuming an excitation temperature of 
20 K. Based on the absence/presence of MIR emission, hence on 
stellar content, we expected a significant spread in relative abun- 
dance between the clumps. Indeed, the N2H^ column density 
coiTelates well with the MIR brightness with the column density 
increasing from clumps 5 to 1 . 

The ratio of C'^0 and N2H^ abundance is visualised in 
Fig. |6] To do this, we have taken the CO and N2H^ abundance 
uncertainties into account. The statistical uncertainties based on 
fits to the line profile is sufficient for N2H^. However for CO, 
while the errors on the fit to line profile is only a few percent, 
the unknown CO excitation temperature, makes it the most 
significant parameter for estimating the uncertainties. We there- 
fore vary the excitation temperature between 15 and 25 K and 
refer to Table |6]for extremes in C'^O column densities. For the 
densities observed in the cloud (> 10^ cm^), Te^^co must be very 
close to the gas temperature. Based on ammonia (1,1) and (2,2) 
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Fig. 6. MOPRA N2H^ moment zero map in contours overlaid on the LABOCA 870;um emission. The red rectangle shows the 
approximate area mapped with MOPRA. Clump labelling is the same as adopted through out the text. The contours have been 
chosen to start at 3 cr, 5 cr in steps of 3 cr where cr - 1 .4 Jy kms"' . 



observations with the Parkes telescope towards clumps 1, 2, and 
3, Wienen et al. (in prep.) find a temperature of <20 K. For 
clumps 2 and 5, we thus adopt a reasonable temperature range of 
10 K to 25 K to calculate the column densities. The mean of the 
largest difference in C'^O column densities at these very differ- 
ent temperatures (Table |6]l is then considered as a conservative 
upper limit for the eiTor on the ratio. A higher temperature ex- 
treme of 35 K is adopted for clump 1 which hosts the hot core 
associated with IRAS 17233-3606. 

Uncertainties in the beam filling factor of N2H^ and C'^O is 
an additional source of error. For the C'^O abundance, we ex- 
pect a value close to 1 for the C'^O column density (because of 
the low critical density of the line), and also for the H2 column 
density given the size estimates obtained in Sect. 14.11 however, 
since we are limited by our single pointing measurements for CO 
and poor resolution for N2H^, the respective beam filling factor 
remains unknown. For the N2H^ abundance, dust emission and 
N2H^ have very similar distributions. Also, the significant differ- 
ences in source sizes extracted from our higher resolution dust 
observations (see TableO indicate that beam dilution might play 
a role. The abundance ratio decreases from 437/214 (clump 1) to 
93 for clump 5 from the brightest end of the filament to the other 



end; however, based on the uncertainties discussed above, the 
observed differences in abundance ratio cannot be ascertained. 
The filament has an average N2H^ abundance (over the three 
clumps) of 2.0x10 '° relative to H2. The derived column density 
and abundanc e is consistent with those found in other infrared 
dark clouds (iRagan et al.ll2006t ISakai et al.ll2008h . 



The fairly constant CO abundance over clumps with very dis- 
tinct levels of star formation suggests that processes other than 
desorption of CO might affect the N2H^ chemistry. In a CS and 
N2H^ study of a small sample of dense cores in high-mass star- 
forming regions, Pirogov et al. (2007) briefly explore alternate 
pathways for the dissociative recombination of N2 H^ other than 
the st andard channel with N2 as the end product dGeppert et al.l 
i2004l) . Howev er, these pathways m ay not be dominant, as previ- 
ously thought (lAdams e t al.'2009'). and new chemical models are 
required to explain this. Observationally, HCO^/H'^CO^ distri- 
bution and abundance might help us in fully ruling out the role 
of CO desorption. 
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Table 8. N2H+ 1-0 line parameters. 



Clump 


vlsr 


Av 






A'n,h+ 










[kms-'] 


[km s-'] 




[K] 


[10'3 cm-2] 


[lO-'O] 


[10-**] 




1 


-3.7 (0.04) 


4.55 (0.06) 


0.2- 1 


11.7-37.3 


5.5 - 11.3 


0.83- 1.68 


3.6 


214-437 


2 


-2.9 (0.01) 


3.06 (0.05) 


1.05 (0.02) 


9.2 (0.4) 


2.5 (0.2) 


2.3 


7.1 


307 


5 


-2.2 (0.04) 


3.36 (0.09) 


1.41 (0.40) 


6.1(1.1) 


1.8(0.8) 


3.7 


3.4 


93 



Note T,o, is the sum of the peak optical depth of the seven hyperfine components, xn^u* Ctc"o) is the N2H^ (C' 'O) abundance. 



6.2. N2H^ line width and velocity field 

We have analysed the line width and velocity fields from the 
NaH^ 3-D data, and the results are shown in Fig. [8] For these, we 
have considered only pixels corresponding to NaH^ integrated 
intensity > 3cr. 

Four clumps are identified in the first moment map (see 
Fig. [8] left panel). These have small differences in velocities 
compared to their line widths (Fig. [8] right panel). The velocity 
distribution for clump 1 is likely complicated by the intense star 
formation activity. If the filament were to be part of a shell-like 
structure (or if it were rotating), one would expect a clear linear 
velocity gradient from one end of the filament to the other Most 
of the extended gas appears to be at a velocity of -3.4 kms 
which we take to be the velocity of the ambient medium. Clumps 
2 and 5 are separated by this lower density extended medium, 
contrary to a linear velocity gradient (see Sect. |7]i. A projection 
effect where clumps 2 and 5 are coincidentally aligned along a 
ring also seems unlikely given that there is no obvious kink in 
the velocity field. 

Of particular interest are two striking features in the veloc- 
ity map; (i) a large-scale (a; 1 parsec) NE-SW velocity gra- 
dient of more than 1.6 kms ' associated with the hot core in 
IRAS 17233-3606 (clump 1), and (ii) a smaller scale velocity 
gradient (» 0.6 parsec) in the northern clump (clumps 5, 6). Both 
velocity gradients are perpendicular to the filament, but in dif- 
ferent directions. Interestingly, the NE(blue)-SW(red) gradient 
towards the hot core coii'el ates well with the directio n and kine- 
matics of the CO outflows ( Leurini et all l2008l 12009 ). Although 
NaH^ is seldom associated with outflows, the interactions of 
outflow shocks with the molecular envelope (traced by NiH^) 
could explain such a gradient towards the hot core (clump 1) 
as detected in the same transition towards the class object 
IRAM 04 1 9 1 4- 1 522 bv iLee et al J (120051) . The axis of the outflow 
(shown in Fig. [8]) would favour this shock scenario. This would 
also explain the very broad lines 5 kms ') observed towards 
the hot core. Such broadening away from the mm core has also 
been observed in lower exc itation NH3 lines in low- and high- 
mass star-forming r e gions (iBachiller et all 119931 : IZhang et all 
II999t iBeuflieT et al.L l2007h . A large-scale infall motion would 
also account for the observed velocity gradient. This large-scale 
infall is o ften observed in regions of high-mass star formation 
dPurcell e t al., 2006). 

The velocities at the peak positions between the northern and 
southern clumps differ by at least > 1.6 kms This is compa- 
rable to the average velocity dispersion (see Fig. [8] right panel, 
and Table [8]l of the clumps within the filament. At this rate, the 
clumps would separate out on a scale of 1.6 pc in 1 Myrs. 

The line width distribution shown in Fig.[8]might reflect dif- 
ferences in the evolutionary status of the clumps along the fila- 
ment. As expected, the hot core (clump 1) exhibits the broadest 
line, and the rest of the filament (clumps 2, 5, and the rest of 
the cloud) has roughly similar line widths with a slightly higher 



peak towards clump 5. Clump 5 is associated with a stronger 
8/im emission than clump 2; however, in all cases the line width 
is highly supersonic. 

Given the line width and the extent of the cloud in projection, 
we can derive the crossing time. Thus the observed turbulents 
motions within the arc of filament would provide an upper limit 
to its lifetime of less than 1 Myr. 

Unlike CO and its isotopologues studied here, the N2H^ tran- 
sition has a high critical density (> 10^ cm-'). Therefore, in 
spite of the lower resolution of the N2H^ data, we are selectively 
probing the dense gas, and not the envelope; i.e., line width-size 
relation does not necessarily apply. This might have the impor- 
tant implication that, although local star formation might be the 
cause of the exceptional line broadening in clumps 1 and 5, an 
external agent might be the cause of the line broadening along 
the whole filament. Interestingly, evidence of large-scale shocks, 
maybe remnants of the IRDC formation process, were found by 
Jimenez-Serra et al. (2010) towards IRDC G035.39-00.33. The 
present data set is insufficient to properly address this scenario. 
Future higher angular resolution observations capable of resolv- 
ing the dense cores structures will allow us to test this scenario. 



7. Discussion 



Recently, iRathborne et al.l (1201 Ol) have studied a sample of 38 
IRDCs in continuum emission at submillimetre and mid- to far- 
IR wavelengths. They found that only a fifth of their clumps 
have 8;um emission, while almost half of them have a 24yum 
point source. In the case of G35 1 .77, the majority of identified 
dust clumps are associated with emission at 24jum, suggesting 
that star formation is going on along the full extent of the fil- 
ament. This is also confirmed by H2 observations at l.lljim 
which detected weak outflows along the full extent of the fila- 
ment (Stanke et al. in prep.). The higher detection rate of 24yum 
co mpact sources in G351 .77 compared to the sample studied 
by Rathborne et al.l (12010') could be due to higher sensitivity 
to low-luminosity objects, given the closer distance of G351.77 
(D < I k pc) compared to the ty pical distance (~4 kpc) of the 
IRDCs of lRafliborne et al.1 (120101) . 

The velocity structure detected in NaH^ shows very strong 
gradients towards individual clumps. This might indicate a gen- 
eral infall within the individual clumps (which might harbour a 
protocluster) or shocks in molecular envelopes. The line widths 
are highly supersonic and might also indicate the presence of 
such shocks. 

The 8/im map (Fig. [3]) shows emission along a semi-shell 
like structure, which is reminiscent of bubbl es of gas and dust 
detected around expanding Hn reg ions (e.g.. lElmegreeiil 119981; 
Deharveng et all l2003Ll200ll2010 '). These structures are poten- 
tial candidates of star formation triggered by massive stars. The 
continuum emission at 870//m is also shaped along a semicir- 
cular structure. Therefore a tempting explanation for the broad 
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Fig. 8. Left Panel: MOPRA NaH^ velocity distribution with dust continuum LABOCA 870/im emission in contours. Red and blue 
arrows indicate the rough direction of the red and blue CO outflow lobes detected towards IRAS 17233-3606 respectively. Right 
Panel: MOPRA N2H^ line width distribution (FWHM) with dust continuum LABOCA 870/im emission in contours. In both maps, 
background image units are in kms and only those pixels that coiTespond to N2H^ integrated intensity > 2.4Jy kms ' have been 
considered significant. 



line profiles associated with dust clumps in G351.77 is shock 
compression due to expansion of Hii regions on pre-existing 
gas. This could also be the triggering event for the intense star 
formation activity shown by several molecular outflows along 
the whole G351.77 (Stanke et al. in prep). We searched the 
SIMBAE0 database for all astronomical objects contained in 
a circle of 0.5° diameter with clump 7 (in the middle of the 
dark filament) at its centre, but could not identify any obvious 
source of trigger. In particular, the only cm cont inuum sources 
in th e region found in the NVSS catalogue dCondon et al.L 
[1998 )) are associated with IRAS 17233-3606 and IRAS 17221- 
3619. From the velocity field derived from the CO observa- 
tions, we know that the LABOCA semicircular structure is as- 
sociated with molecular gas at two different velocities, one at ~ 
-3 km (IRAS 17233-3606, G351.77 and clump 3), the other 
at ~ -22 km s"' (IRAS 17221-3619, IRAS 17227-3619 and the 
position G351.62). The near kinematic distance of IRAS 17221- 
3619 and IRAS 17227-3619 is ~ 3.6 kpc(DFAR ~ 13 kpc), while 
the two velocity components of the position G35 1 .62 correspond 
to kinematic distances of ~ 600 pc and -3.6 kpc (Z)far ~ 16 
and 13 kpc, respectively), derived by using the r otation curve 
ofjBrand & Blitz (1993). For IRAS 17227-3619. iPeeters et all 
( I2OO2I) derived a kinemat ic distance of 3.4 kpc based on IR re- 
combination lines, while iBorissova et al estimated dis- 
tances between 1.9 and 3.2 kpc based on H- and K-band spectra 
of the source. Thus, the semicircular structure detected at 870/im 
in the LABOCA map is likely a projection effect of two different 
molecular clouds. However, from our pointed CO observations 
we can neither infer the true distribution of the two cloud com- 
plexes nor verify whether the western half circle, delineated by 
diffuse 8//m emission (see Sect. |3]l, is connected to one of these 
two clouds. 



8. Conclusions 

In summary, we have characterised the G351.77 dense filament 
in thermal dust continuum, NaH^ (1-0) and CO (various iso- 
topologues at low excitation) emission. The molecular line emis- 
sion traces the thermal dust continuum very well. Three main 
clumps in N2H+ and eight clumps (at higher resolution) in the 
dust continuum were observed. These clumps also have asso- 
ciated emission from CO isotopologues. We analysed the inte- 
grated intensity, velocity, and line width distribution along the 
filament. The column density and mass of the clumps were 
also derived. The velocity structure shows very strong gradi- 
ents across two clumps (1 and 5). This might indicate either a 
global infall within the individual clumps (which might harbour 
a protocluster) or shocks in molecular envelopes due to outflow- 
envelope interaction. All the clumps have supersonic line widths. 
Such supersonic line widths along the filament might also indi- 
cate the presence of shocks. Clearly, the most evolved clump 
is the hot core that shows (i) the brightest MIR emission, (ii) 
broadest line width, (iii) largest velocity gradient, and the (iv) 
most massive and dense core. Most mm clumps have a mass 
> 10 Mq. 

Our study shows that G35 1 .77 has ongoing star formation in 
different evolutionary stages along the filament. Therefore, given 
the relatively close distance of this object compared to typical 
sites of massive star formation, G351.77 is an ideal target for 
investigating the role of IRDCs in the early phases of the for- 
mation of massive stars and stellar clusters, and for studying the 
population of low-mass stars in these objects. 
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